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(57) ABSTRACT 

A computer-implemented model of fan section of a gas 
turbine engine accounts for the turbulence in the gas flow 
emanating from the rotor assembly and impinging upon an 
inlet to the stator vane cascade. The model allows for 
user-input variations in the sweep and/or lean angles for the 
stator vanes. The model determines the resulting acoustic 
response of the fan section as a function of the turbulence 
and the lean and/or sweep angles of the vanes. The model 
may be embodied in software that is rapidly executed in a 
computer. This way, an optimum arrangement in terms of 
fan noise reduction is quickly determined for the stator vane 
lean and sweep physical positioning in the fan section of a 
gas turbine engine. 

12 Claims, 9 Drawing Sheets 
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STAGGER ANGLE 9 
Vane sections at constant radius 
unwrapped ontox* y d plane (duct 
coordinates). 

x& y 0 coordinates aligned with vanes 
at stagger angle 6. z d and z a axes 
aligned with leading edge of reference 
vane. 




Other views in the x^ y^ z 0 system. 
For the lean and sweep operations to 
follow, consider the vanes attached to 
the y d axis with “hooks”. Coordinate 
systems follow the reference vane; 
vane #1 stays parallel to ref. vane 
(hooked dty d = g). 


LEAN ANGLE w _ 

Lean the reference vane about the x 0 
axis through angle ifr z'axis follows 
leading edge, x 'axis coincides with x„ 
axis. Because of the “hook”, the gap 
in the y d system is not changed. 


SWEEP ANGLE y/ s 
Sweep the vane in its own plane 
through angle y/ r Velocity vector 
does not change. c d is die chord 
measured at constant radius in the fan. 
c = c d cos y/ z 


The x,y, z coordinates are now 
Glegg’s cascade coordinate system. 
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'Comparison with ADP#1 (722ft /sec) data scaled to 130 inches ' 

.80 35. 0. 0. 45 G theta PsilD PsisD Vanes 

.5006 0.035 0.02 Mr L I 

2223. 4.6 2.203 1125. Po R DelR Ao 

3 Spectrum type: 1 for const BW, 3 for 1/3 OB 

100. 6300. ** See below for this line 

.80 .80 DeltaAlf DeltaNu 

** For const BW, input: fl delfreq Nfreq BW 

** For 1/3 OB, input: first freq, last freq (e.g. 63. 500.) 
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PREDICTING BROADBAND NOISE FROM A 
STATOR VANE OF A GAS TURBINE ENGINE 


GOVERNMENT RIGHTS 

The U.S. Government may have rights in this invention 
pursuant to NASA contract NAS3-27727. 


TECHNICAL FIELD 


This invention relates generally to predicting the amount 
of noise in a gas turbine engine, and more particularly to a 
computer-implemented method that predicts the amount of 
broadband noise from the stator vane cascade in a gas 
turbine engine. 


BACKGROUND ART 
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An aircraft gas turbine engine includes a fan section, a 
compression section, a combustion section, and a turbine 
section. An annular flow path for the working medium gases 
extends axially through the engine. The fan section has a 20 
rotor assembly that includes an array of rotor blades that are 
angled with respect to the approaching gas flow. The fan 
section also has a stator assembly that includes an array or 
cascade of stator vanes. Each vane is typically arranged in 
a radial direction outward from the center axis of the engine. 25 
The stator vane cascade is disposed axially downstream of 
the rotor blade array in the gas flow path. The vanes receive 
and guide the direction of the flow of gases exiting from the 
upstream rotor blade array. The stator vanes connect 
between an outer duct wall and an inner duct wall within the 30 
engine. The duct walls extend circumferentially with respect 
to the flow path to form the boundaries for the working 
medium gases in the fan section. 

As the rotor assembly rotates, the blades do work on the 35 
gases to increase the pressure of the gases. The rotor blades 
also increase the velocity of the gases and direct the flow of 
gases from the engine axial direction to the blade rotation 
direction. The gases are then flowed past the rotor blade 
array to the stator vane cascade, which redirects the flow of 4Q 
gases to the axial direction. By reorienting the flow in this 
manner, the stator vane cascade increases the recovery of the 
flow energy of the gases into thrust. 

As the working medium gases travel along the engine 
flow path, the gases are pressurized in the rotating fan and 45 
compression sections, which causes the gas temperature to 
rise. The hot, pressurized gases are burned with fuel in the 
combustion section to add energy to the gases. The gases are 
then expanded through the rotating turbine section to pro- 
duce useful work for pressurizing the gases in the fan and 50 
compression sections. The expanding gases also produce 
thrust for propelling the aircraft forward. 

The gas flow through the engine generates acoustic 
energy or noise. Aircraft engine manufacturers are con- 
cerned with the adverse effect of excessive noise levels on 55 
passengers, aircraft personnel, and residents in close prox- 
imity to airports. Due to increasingly stringent noise restric- 
tions placed upon aircraft that operate in certain geographic 
areas and at certain times, there is a persistent need for 
quieter aircraft engines. 60 

The principal sources of noise in an aircraft gas turbine 
engine are jet or exhaust noise, core noise, and fan noise. Jet 
noise results from mixing of the high-velocity engine 
exhaust gas stream with the ambient air. A considerable 
amount of turbulence is generated when these two streams, 65 
which are traveling at different velocities, mix together. This 
turbulence generates jet noise. In a turbofan engine, there are 


two exhaust streams; therefore, there are two sources of jet 
noise. One noise source is the turbulent mixing of the fan 
exhaust stream with the ambient air. The other noise source 
is the turbulent mixing of the engine core exhaust stream 
with the fan exhaust stream and the ambient air. 

Core noise consists of compressor noise, combustion 
noise and turbine noise. Compressor and turbine noise are 
caused by the unsteady blade forces and fluid stresses when 
fluids are compressed for driving the turbines. Combustion 
noise results from the turbulence generated by the burning of 
fuel in the combustion chamber. 

Fan noise is often the predominant noise source in a 
high-bypass ratio turbofan engine. Fan noise is caused by 
non-uniform gas flow exiting the rotor blades and impinging 
upon the stator vanes in the fan section. As each rotor blade 
passes through the gases, the blade leaves a wake or track of 
turbulent gases behind it at the trailing edge of the blade. 
This wake is commonly referred to as the rotor wake 
turbulent flow. Also, tracks of non-turbulent flow exist 
between the rotor blade array and stator vane cascade in the 
areas not directly behind the blades. The turbulent flow 
impinges on the stator vanes and generates much of the noise 
radiated by the engine. This fan noise is also commonly 
referred to as rotor/stator interaction noise. 

In addition, secondary flow patterns exist adjacent the tips 
of the rotor blades due to the interaction of the blade tip with 
the boundary layer of the engine outer duct wall. This 
interaction introduces further turbulence (i.e., “tip vortex”) 
into the wake turbulent flow at the blade tip region. Also, 
turbulent flow emanates from the hub or root portion of the 
rotor blade (i.e., “hub vortex”). The wake turbulent flow 
(i.e., “downwash”) from rotor blades sweeping past stator 
vanes produces pressure fluctuations on the vane surfaces. 
These fluctuating aerodynamic pressures on the vane sur- 
faces produce forces that generate noise. 

The rotor blade wake turbulent flow has a steady com- 
ponent and a random component. The steady component is 
also commonly referred to as the harmonic or periodic wake 
component. 

The random component is also commonly referred to as 
the broadband wake component. The random component is 
represented by the turbulent kinetic energy, which varies 
dramatically across the span of the stator vane inlet. The 
turbulent kinetic energy is often greater at the root or hub 
and tip regions of the rotor blade as opposed to the mid-span 
region of the rotor blade. The turbulent kinetic energy 
present in the blade root region is typically attenuated as it 
is absorbed downstream in the low-pressure compressor 
section of the engine. 

The harmonic and broadband components of the wake 
turbulent flow are related to the noise spectrum components. 
The harmonic wake component causes harmonic or tonal 
noise and the broadband wake component causes broadband 
noise. Tonal noise is noise at specific frequencies, which are 
multiples of the rotor blade passage frequency. This tonal 
noise has a distinct sound that can be heard above the 
background noises. The amount of tonal noise generally 
depends upon the number of stator vanes and rotor blades, 
the geometry of the duct walls bounding the flow path for the 
working medium gases, the velocity of the gases, and the 
rotor speed. On the other hand, broadband noise is distrib- 
uted over a wide range of frequencies, rather than being at 
specific, discrete frequencies. The noise signatures of 
modem turbofan engines tend to be dominated by broadband 
noise. 

It is difficult to suppress or attenuate fan noise because of 
the interdependence of the mechanisms that contribute to 
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this noise and the basic aerodynamic operation of the fan 
section of the engine. The prior art contains noise suppres- 
sion structures adapted specifically for retrofit or original fit 
on an aircraft gas turbine engine. Typically, the noise sup- 
pression structure consists of sound attenuating liners 5 
applied to the nose cowl, the nose dome and the fan flow 
path components of the engine. In typical constructions, the 
sound absorption material lines the inlet duct and nozzle of 
a turbojet or turbofan engine to suppress the noise generated 
within the flow path. However, significant aerodynamic 
losses (e.g., thrust) result from the addition of noise sup- 
pression structures that provide for acceptable levels of tonal 
and broadband noise. 

Another known method of attempting to reduce the fan 
noise involves selecting blade/vane ratios to satisfy the cut 15 
off criterion for propagation of noise at the fundamental 
rotor frequency. It is also known to increase the axial 
spacing between the rotor assembly and the stator assembly 
to reduce fan noise. 

The aerospace community is aware of the potential for 2 o 
reducing the tonal noise component of the turbulent rotor 
wake by adjusting the angular physical positioning of the 
stator vanes in one or both of two different angular direc- 
tions. This angular physical positioning is known as “lean” 
and “sweep”. Vane sweep is the axial displacement of the 2 S 
vane such that the tip region of the vane is disposed farther 
downstream axially than the vane hub. Correspondingly, 
vane lean is a circumferential displacement of the vane 
stacking line relative to the radial direction of the vane. Lean 
of the vane tips is normally in the direction of blade rotation. 30 
Sweep and lean vane positioning may be performed inde- 
pendent of each other, or they may be done simultaneously. 
Swept and/or leaned stator vanes reduce tonal noise by 
reducing the severity of the rotor wake interaction with the 
vanes. The axial spacing (i.e., upstream, downstream) of the 35 
rotor assembly with respect to the stator assembly may also 
typically be varied with various sweep and/or lean configu- 
rations. 

Recent empirical laboratory testing of an engine fan 
section has revealed the unexpected potential for use of 40 
stator vane sweep and/or lean physical positioning to reduce 
the broadband noise component of the rotor turbulent wake 
flow. See Woodward, R. R, et al., “Benefits of Swept and 
Leaned Stators for Fan Noise Reduction”, AIAA Paper 
99-0479, presented at the 31 th Aerospace Sciences Meeting 45 
and Exhibit, 1999, Reno, Nev., Jan. 11-14. Therein it was 
described that sweep only, and sweep together with lean, 
were both effective for reducing the amount of broadband 
noise at the stator vane inlet. The quantity of measured 
broadband noise was given as a function of various lean and 50 
sweep angles. The acoustic response results revealed sig- 
nificant reductions in both tonal noise and broadband noise 
beyond what could be achieved through the conventional 
approach of increasing the axial spacing between the rotor 
assembly and stator assembly. 55 

Yet, instead of attempting to determine optimum vane 
sweep and/or lean positioning through the trial and error 
empirical testing approach, it is desired to accurately model 
various vane sweep and lean angles and the resulting acous- 
tic response of the fan section. This is done in the present 60 
invention to accurately predict the amount of broadband 
noise for various sweep and/or lean angles. The model may 
be embodied in software that is rapidly executed in a 
computer implementation of the model. This way, an opti- 
mum arrangement in terms of fan noise reduction is quickly 65 
determined for the stator vane lean and sweep physical 
positioning in the fan section of a gas turbine engine. 
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DISCLOSURE OF INVENTION 

An object of the present invention is to provide a 
computer-implemented method for accurately predicting the 
amount of broadband noise at the inlet of a stator vane 
cascade of a gas turbine engine. 

According to the present invention, a model of a fan 
section of a gas turbine engine accounts for turbulence in the 
gas flow emanating from the rotor assembly and impinging 
upon the stator vane cascade at the inlet thereof. The model 
also allows for variations in the sweep and lean angles of the 
stator vanes. The model determines the resulting acoustic 
response of the fan section as a function of the turbulence 
and the lean and/or sweep angles of the vanes. The model 
may be embodied in software that is rapidly executed in a 
computer or workstation. The user enters desired values for 
the turbulence and the vane lean and/or sweep angles. The 
computer executes the program instructions embodying the 
model to solve the equations for the acoustic response of the 
stator vane. This way, an optimum arrangement in terms of 
fan noise reduction is rapidly determined for the stator vane 
lean and sweep physical positioning in the fan section of a 
gas turbine engine. 

The present invention has a practical application and 
utility in that it provides for a computer-implemented tool to 
assist gas turbine engine designers in designing the fan 
section. The present invention allows the designer to “vir- 
tually design” certain aspects of the stator vane cascade 
using a computer model, to achieve a minimum or optimum 
amount of fan noise emanating from the engine. Specifically, 
it allows the designer to vary the lean and/or sweep angles 
of the stator vanes, and then to assess the effect of those 
angles on the amount of fan noise generated by the inter- 
action of the rotor wake turbulence with the stator vane 
cascade. This type of “what if’ analysis can be performed 
relatively much more accurately, rapidly and inexpensively 
in a virtual computing environment, as compared to a 
traditional laboratory empirical “trial and error” approach. 

The above and other objects, advantages and features of 
the present invention will become more readily apparent 
when the following description of exemplary embodiments 
of the present invention is read in conjunction with the 
accompanying drawings. 

BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 is a perspective view, partially cut away, of a gas 
turbine engine of the axial flow, turbofan type having a rotor 
blade array and a stator vane cascade; 

FIG. 2 is a side view of a rotor blade and a stator vane 
within the engine of FIG. 1, and illustrating the various 
components of the turbulent gas flow in the rotor blade wake 
as it impinges upon the stator vane; 

FIG. 3 is a block diagram of a computer workstation 
having a central processing unit for executing instructions 
that implement the noise prediction model of the present 
invention; 

FIG. 4 includes various graphs and a diagram that illus- 
trate various geometrical components of the gas flow of FIG. 
2 with respect to a stator vane cascade of FIG. 1; 

FIG. 5 includes various graphs that illustrate the sequence 
of rotations that relate cascade coordinates to duct coordi- 
nates with respect to stator vanes in the cascade of FIG. 1; 

FIG. 6 includes graphs that illustrate exemplary lean and 
sweep angles for the stator vanes in the cascade of FIG. 1; 

FIG. 7 includes graphs that illustrate lean and sweep 
angles as defined via front and side views of stator vanes in 
the cascade of FIG. 1; 
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FIG. 8 include graphs that illustrate the turbulent gas flow 
from the rotor blade array impinging on stator vanes in the 
cascade of FIG. 1; and 

FIG. 9 is a listing of a computer input file containing 
exemplary values for various parameters used by a computer 
program executing on the workstation of FIG. 3 in deter- 
mining the acoustic response of the stator vane cascade of 
FIG. 1. 

EXEMPLARY EMBODIMENTS OF THE 
PRESENT INVENTION 

FIG. 1 illustrates a gas turbine engine 10 of the axial flow, 
turbofan type having an axis of rotation, A r . The engine 10 
may comprise the model PW4000 commercial aircraft gas 
turbine engine, provided by the Pratt & Whitney division of 
United Technologies Corporation, the assignee of the 
present invention. A nacelle 12 surrounds the engine 10 and 
is adapted to both support and position the engine from a 
structure such as an aircraft wing (not shown). The engine 10 
includes a fan section 14, a compressor section 16, a 
combustion section 18, and a turbine section 20. A primary 
annular flow path 22 for working medium gases extends 
through these sections 14—20. A secondary annular flow path 
24 for working medium gases is disposed radially outward 
of the primary flow path 22. 

The engine 10 includes a stator assembly 26 extending 
circumferentially about the engine axis of rotation. The 
engine 10 also includes a rotor assembly 28 disposed about 
the engine axis of rotation. The rotor assembly 28 includes 
a rotor disk 30 and an array of rotor blades 32. Each rotor 
blade 32 has a root or hub region 34, a mid-span region 36, 
and a tip region 38. Each rotor blade 32 extends outwardly 
from the rotor disk 30 across the working medium flow paths 
22, 24. 

The engine 10 also includes an outer duct wall 40 that 
extends circumferentially about the engine axis of rotation. 
The outer wall 40 outwardly bounds the secondary annular 
flow path 24 for the working medium gases. The tip region 
38 of each rotor blade 32 is in close proximity with the outer 
wall 40. An inner wall 42 (FIG. 2) is spaced radially 
inwardly from the outer wall 40. The inner wall 42 extends 
circumferentially about the engine axis of rotation and 
inwardly bounds the secondary annular flow path 24 for the 
working medium gases. 

The stator assembly 26 includes an array or cascade of 
stator vanes 48 spaced axially downstream from the rotor 
blade array 28. Each stator vane 48 extends radially between 
the outer and inner duct walls 40, 42 and across the 
secondary annular flow path 24. A stator vane 48 has a radial 
span that may be measured at the leading edge of the vane. 

FIG. 2 is a simplified side view of a portion of the 
turbofan engine 10 illustrating a rotor blade 32 and a stator 
vane 48 disposed in the secondary annular flow path 24. The 
stator vane 48 is illustrated in a common radial orientation 
connected between the outer and inner duct walls 40, 42. 
Also illustrated are the elements of the rotor blade turbulent 
wake, including the tip vortex or wakes 50, the med-span 
wakes 52, and the hub or root vortex 54. These non-uniform 
turbulent flow components in the rotor blade wake impinge 
on the stator vanes 48, and may generate most of the noise 
radiated by the engine 10. 

The wake includes a periodic component, which gener- 
ates harmonic noise, and a random component, which gen- 
erates broadband noise. The present invention focuses 
exclusively on predicting the broadband noise component 
for various angular positions of the stator vanes 48. The 
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random flow is considered to be turbulent, ignoring any 
random instability effects. Generally, the rotor exit flow is 
highly inhomogeneous because the turbulent blade wakes 
50-54 are not mixed and fully merged with each other. Also, 
5 bands of turbulence exist at both the outer and inner walls 
40, 42. The statistics of this turbulence, like those of the 
mid-span wakes 52, exhibit a blade- to -blade periodicity. 

An analytical model is developed herein for predicting the 
amount of noise generated by the rotor wake turbulence at 
10 the inlet to the stator vane cascade. The model treats rotor 
wake turbulence via three-dimensional spectra and covari- 
ance functions. Special attention is paid to effects of lean and 
sweep angular positioning of the stator vanes and to inho- 
mogeneity of the turbulence. The result of this rigorous 
15 analysis is several tractable mathematical equations that 
predict the sound power spectra upstream and downstream 
of the stator vane cascade as an integral over the wavenum- 
bers of the incoming rotor wake turbulence field. The sound 
power equations are also functions of various numerical 
20 values for stator vane sweep and/or lean angles. The equa- 
tions are program coded in software for execution on a 
computer for solution thereby. In an exemplary embodiment, 
the equations are coded in Fortran and run on a commer- 
cially available workstations or personal computer 70. 

25 Referring to FIG. 3, there illustrated is a simplified block 
diagram of the workstation 70 available from, e.g., Sun 
Microsystems or Silicon Graphics. The workstation 70 
includes a central processing unit (“CPU”) 72, which may 
comprise the Intel Pentium microprocessor. The workstation 
30 70 also includes various types of memory 74 for storing 
different types of information and data, including, for 
example, instructions executed by the CPU 72 and the 
results of calculations performed by the CPU. The various 
known types of memory 74 may comprise read only 
35 memory (“ROM”), random access memory (“RAM”), and 
electrically programmable read only memory (“EPROM”). 

The workstation 70 also includes the means 76 for com- 
municating with external storage devices 78. These devices 
40 78 typically include various known types of storage media, 
including floppy disks and CD-ROMs. The devices 78 may 
store the underlying functional instructions and accompa- 
nying data that comprise the noise prediction modeling 
software of the present invention. The stored instructions 
45 may be executed by the CPU 72. The means for communi- 
cating 76 with these devices 78 handles the reading and/or 
writing of information to and from these external storage 
devices 78. 

A visual display device 80, such as a CRT monitor, is 
50 typically included as part of the workstation 70. The monitor 
80 displays various types of visual information to the 
workstation user, including graphical data illustrating the 
amount of noise at the inlet of the stator vane cascade for 
various lean and/or sweep vane angles. The CPU 72, 
55 memory 74, external storage communication means 76, and 
visual display device 80 may communicate with each other 
via a bus 82, which comprises a plurality of communication 
signal lines. 

In operation of an exemplary embodiment of the present 
60 invention, the analytical model described in detail herein- 
after yields one or more equations that calculate the amount 
of noise at the inlet of the stator vane cascade. The equations 
are coded in the Fortran software language, although any 
suitable programming language may be used. The user then 
65 inputs into the computer program values for various param- 
eters (e.g., turbulence, stator vane lean and/or sweep angles) 
needed to calculate the acoustic response of the cascade. The 
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CPU 72 performs the calculations, and a visual output of the 
results of the calculations is provided on the display screen 
80. The use may then modify one or more of the various 
parameter values and re-execute the calculations. This can 
be done until an optimum or minimum desired value for the 
amount of stator vane inlet noise is achieved. 

With respect to the analytical model developed herein, the 
only simplification regarding representation of the rotor 
wake turbulence is use of Taylor’s hypothesis (i.e., the 
frozen gust assumption). This reduces the turbulence repre- 
sentation from a four-dimensional wavenumber-frequency 
spectrum to a three-dimensional wavenumber spectrum and 
results in manageable noise equations. The principal sim- 
plification regarding geometry is use of a rectilinear stator 
vane cascade model. At low frequencies, this may be prob- 
lematical. However, at mid and high frequencies where 
noise generation is dominated by small scale eddies, the 
rectilinear cascade approximation should be acceptable. 

The relatively complex analytical derivation of the noise 
prediction model of the present invention starts with a 
known, wake impingement analysis model. See Glegg, S. A. 
L., “The Response of a Blade Row to a Three Dimensional 
Gus”, Florida Atlantic University Report, September 1996. 
This model is hereafter referred to as the “Glegg model”. 
The subsequent mathematical analysis presented herein 
develops the sound power equations from the Glegg model. 

The Glegg model essentially comprises a flat plate, rec- 
tilinear cascade acoustic response theory. This theory 
includes a spanwise flow component for treatment of the 
sweep effect, and a spanwise wavenumber component, in 
addition to the usual streamwise and gapwise wavenumber 
components. Glegg’ s theory is for waves that are sinusoidal 
in space and time. The following mathematical analysis 
connects Glegg’s theory to the turbulence interaction prob- 
lem. The resulting equations give the spectra of upstream 
and downstream sound power generated by rotor exit wake 
turbulence. Inhomogeneity is dealt with in terms of integrals 
of the turbulence covariance function, or of the turbulence 
spectrum function, at the stator cascade face. Use of the 
resulting noise prediction model theory of the present inven- 
tion is described and illustrated herein by a simple applica- 
tion of the Liepmann turbulence spectrum for homogeneous 
turbulence, in place of the general turbulence spectrum. This 
teaches one of ordinary skill in the art how to integrate over 
a turbulence field and how to transform wavenumbers for 
turbulence represented in different coordinate systems (i.e., 
the “cascade” and duct coordinate systems described in 
detail hereinafter). 

For background information on the Glegg model, refer- 
ence is made to FIG. 4, which illustrates four stator vanes 48 
in the cascade (FIG. 1). Geometry is constant in the z 
direction and the background flow is uniform, U=(U, 0, R). 
The stator vanes 48 are shown as unloaded flat plates. 
Cascade gap (i.e., the distance between vanes 48), chord, 
and stagger angle (with respect to the y axis) are indicated 
thereon as s, c, and X, respectively. The unsteady flow is 
harmonic in space and time with upwash given by: 

w(x, t)=w 0 e i ^ y ° x+CLy+vz ~ <1 ^ (1) 

This equation represents a plane wave that is harmonic in 
time with frequency co/2jt and upwash complex amplitude 
w 0 . It is also harmonic in space with x, y, and z wave 
numbers equal to y 0 , a, and v. By use of Wiener-Hopf 
analysis, Glegg derived an equation equivalent to the fol- 
lowing for the velocity potential of the scattered acoustic 
waves. 


d ± (x t) = + 7rw ° c2 V 1 & D (K) i [-A± (x-yd/h) + (<r-2xk )y/h+vz] -tut 

5 

Where 



M = U/a jS = V 1 - M 2 (x) g =a)-Wv 

(3) 

10 

s e = Vd 2 + (Ph 2 tan^ e = dj fih Xf = kM + rjf 

(4) 


£ = - (nff t = (v IP) 1 k = u> g Kafi) 

(5) 

15 

nt = ~fk sin^ e ± cos^ ^k 2 - f k =(o--2xk+ xMd) / s e 

(6) 


and a=y 0 d+ah is the interblade phase angle. D is the Fourier 
transform of the discontinuity in potential across the blade 
32 and wakes (in the form of an infinite product), and is the 
primary result of Glegg’s derivation. The velocity, pressure, 
and density perturbations associated with the acoustic wave 
can be obtained from Equation 2 via: 

w= V §p=- poDty/Dtp '=p! a 2 (7) 

25 

where p 0 and a are the ambient density and speed of sound, 
respectively. Glegg’s notation is adhered to closely. 
However, Glegg used of for the radian frequency (herein it 
is to), co for the shifted frequency (herein it is cog), and m for 
30 the scattering index (herein it is k). Also, the D function of 
this report is non-dimensional; thus, to obtain the D function 
of Glegg’s report, multiply the non-dimensional version by 
w 0 c 2 . 

The equations above yield the acoustic waves scattered by 
35 the stator vane cascade for a single, planar wave input. 
Scattering index k runs over an infinite range, but, as usual 
in this kind of formulation, only a finite number of waves are 
cut on (propagate undiminished); the remaining waves 
decay exponentially and, thus, carry no acoustic energy. 
40 Cuton is governed by the argument of the square root 


yj^e-fk ■ 

45 When the frequency is high enough, the argument is positive 
and the waves are cut on. The exponential dependence on the 
space and time variables in Equation 1 permits treatment of 
any inflow field via standard Fourier transform methods, as 
shown as follows. 

50 The analysis in support of the development of the noise 
prediction model of the present invention requires four 
coordinate systems: a first is aligned with the fan duct for the 
lean/sweep definition; a second is aligned with the mean 
flow and the fan radius for the generalized turbulence 
55 description; a third is aligned with the mean flow and the 
vane normals to apply a classic upwash turbulence spec- 
trum; and a fourth is aligned with the cascade geometry to 
apply Glegg’s theory. The transformations and lean/sweep 
conventions are defined herein, and the most important 
60 geometrical relations between the chord, gap, and stagger in 
the duct coordinate system and the cascade coordinate 
system are summarized. Lean and sweep are also defined in 
two different coordinate conventions. 

The coordinate systems are shown in FIG. 5. Stagger, 
65 lean, and sweep angles are defined via an ordered sequence 
of rotations about coordinate axes as follows. The top of 
FIG. 5 represents the stator vane 48 cut at constant radius 
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and unwrapped onto the plane of the paper. The z d and z 0 
axes are coincident with each other and with a fan radius. 
Mean flow U 0 is aligned with the x 0 direction and has no 
radial component. In the transformations that follow, the 
reference vane (vane #0 which passes through the origin) 5 
and vane #1 are handled differently. To visualize the 
rotations, the reference vane is rotated and then the coordi- 
nate system is rotated to re-align the z axis with the leading 
edge of the reference vane. 

The other vanes 48 have hooks on their leading edges 10 
attaching them to the y d axis with a separation g that does not 
change. Thus, throughout the lean and sweep rotations, the 
vanes 48 always penetrate the constant radius plane (z o= 0) at 
the location and orientation given by g and 0 at the top right 
of FIG. 4. 15 

Referring also to FIGS. 6 and 7, lean is an out-of-plane 
rotation of the vane 48 about the x 0 axis through angle 
Lean does not change the component of inflow mean veloc- 
ity normal to the leading edge. Sweep is an in plane rotation 
of the vane 48 about the y' axis through angle ty s . Sweep 20 
reduces the flow component normal to the leading edge by 
the factor cos ty s . Relationships between the velocity com- 
ponents can be deduced from FIG. 5, as can the relation 
between the chords c and c d . 

Glegg’s stagger dimensions in FIG. 4 are given by: 25 

^=(sin 0 cos ip s -cos 0 sin ip/ sin ip s )g 

h={ cos 0 cos i^)# (8) 

From these, Glegg’s gap and stagger angle are given by: 3Q 

s = V d 2 + h 2 (9) 



The chord relation is: (11) 

Where C d is chord measured atconstant radius (along the 
flow direction). The velocity components are related to the 
total mean inflow velocity U 0 by 40 


U=U 0 cos 

W=U 0 sin (12) 

The lean and sweep angles can also be defined by a 
convention more suitable to viewing on design drawings. In 45 
FIG. 7, the stator vane cascade is viewed along the fan axis 
of rotation for lean and from the side (similar to FIG. 2) for 
sweep. 

The lean and sweep angles have “hats” in the following 
equations to distinguish them from the corresponding lean 50 
and sweep angles in FIG. 5. The relationship between the 
two systems is given by: 

tan ijj^cos 0 tan ij>/-sin 0 tan 

tan i|j s =cos ^/(sin 0 tan ij^+cos 0 tan (13) ^ 

Thus, the lean and sweep angles in the “duct system” in FIG. 

7 can be converted to the “no hat” system for use with 
Equations 8-12. Equations 13 must be solved sequentially. 

Glegg’s theory was written for three-dimensional, planar, 60 
harmonic waves. Glegg’s theory can be extended to an input 
having any waveform and spatial distribution via the fol- 
lowing Fourier transform: 

w(x, t)=\\ W(K,<ti)e i(K - x - ,JJt )dKd<ti (1 4) 

65 

where here, and throughout this mathematical analysis, 
integration limits are from -ooto oo) unless otherwise speci- 


10 

fled. Vector wavenumber K is shorthand for (y 0 ,a,v). For 
discussion hereinafter, the inverse of Equation 14 is: 

W(K, (o) = w(x, t)e- i(K ' r ~ ajt) dxdt 


To insure convergence of this integral, gust velocity w is 
considered to be non-zero for a large, but finite, block of 
fluid and for the associated time, T, required for that flow to 
pass the stator vane cascade at speed U 0 . Later, T and the 
fluid volume will approach infinity as a limit. The limiting 
process will become more precise as the derivation pro- 
ceeds. With the application of Equation 14 to Equation 2, 
Glegg’s potential can be generalized to: 


(j)=(^ t)=J <3>=(x,G))e ^‘d co 


(16) 


where 


^(X, CO) = 


( 17 ) 


± 



£ i[~ ( x-yd/h)+(<r-2Kk )y/h+vz] ^ g 


It is desired to compute the spectrum of sound power 
scattered by the stator vane cascade. The starting point is the 
expression for the acoustic energy flux vector applicable 
It is desired to compute the spectrum of sound power 
scattered by the stator vane cascade. The starting point is the 
expression for the acoustic energy flux vector applicable to 
waves in a uniformly moving medium with mean properties 
given by the density p 0 , speed of sound a, and velocity U 
(see Goldstein, M. E., “Aero acoustics”, McGraw-Hill, New 
York, 1976). 


/ = | — + U -wj(pow+ Up') 


( 18 ) 


This is the time dependent power per unit area in terms of 
the acoustic pressure, density, and velocity, p, p', and u, 
respectively. The desired form for spectrum of sound power 
flux based on velocity potential is: 


-ilnpo r u 

I 0J = — w<t>* V4>+ -^-(«w<t> - U ■ 

T L a 2 


(19) 


This is the local (dependent on position variables) power 
flux vector per unit frequency. The following equation is 
used to find the acoustic power leaving the stator vane 
cascade (i.e., downstream of the cascade): 

I^=I^fr (20) 

Then, integrate over the cascade frontal area. fr=±(-h/s, d/s, 
0 ) are the unit normal vectors for a plane parallel to the 
cascade leading edge. The signs +/- correspond to normals 
on the upstream/downstream sides of the plane. By applying 
Equation 19 to Equation 17 and performing the indicated dot 
product, the result is: 

± -ihtpo 71 2 c 4 cl) ( 21 ) 

^ = T We 
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W'(K, io)W(K', to) 


-continued 

ke -ft V < 2 “ /3 


X 



.S'/I 


(<x / - 27TA:') - 


h U , + 

__ (w + ty A p-vW)j|x 

c ‘K A A -A £/ H(c-27i* Mo J ' )]}#-( V-v' )z} rf 


W(K , w) = p~yff w ( x ’ t)e~ i{K - x -° Jt) dxdt 


(24) 


5 

wherein the wave number and space variables are expressed 
explicitly. In the x 0 , y 0 , z 0 system, the wave number vector 
is written as K=(k^, k y , k z ) so that K*x becomes k^.x 0 +k y y 0 + 
k^z 0 . The velocity field in the fluid-fixed system is written as 
10 w(x, y 0 , z 0 ) where, by Taylor’s hypothesis, it is time 
independent. With this notation, w can be written: 


where now the summations are over the limited ranges of k 
and k 1 corresponding to cut on (energy carrying) modes. 
From Glegg’s analysis, it can be shown that the term in curly 
brackets { } is equal to 


w(x y t)=w(x 0 - U 0 t,y 0 ,z 0 ) 


(25) 


15 


Substituting into Equation 24 and changing integration 


variables via x=x 0 -U 0 t lead to: 


Hence, Equation 21 reduces to: 


W(K, oj) = (26) 

20 JL J J J w(x, y 0 , zo)e- i{kx * +k yy° +kzZo) ^ j ' 

dxdy 0 dzo 


t = - 


fiss e 


-f f— zz 

j j k k> 


QtfD-atmp) (22) 


t t ' V*?-/* 2 

c ‘K A A ~ X 'k> Y x ~y d l h )-[(o--2b* )-( cr / -2 nk ’ )]y//z-(v-v' )z}^ ^ 


25 The expression in square brackets [0] is a delta function 
expressing the result of Taylor’s hypothesis that frequency, 
stream-wise wave number, and convection velocity are 
related through co=k^.U 0 . Thus: 


where S ww , is the upwash spectrum: 


. 2 K . 

S WW (K, K\ oj) = y r(/ <, (o)W{K', oj) 


(23) 


30 W(K, oj) = 

6(d), j' j " J yo , Zo )e~ t{kxX+k y y o + ^ z ° ] dxdy 0 dzo 


(27) 


Equations 22 and 23 are the most general form of the sound 35 
power spectrum. Equation 22 represents power per unit 
frequency per unit area of stator vane cascade face. It has 
been averaged over time but is still a function of the space 
variables x, y, z. However, the equations are generally too 
unwieldy for numerical work, considering the six-fold inte- 40 
gral and double sum. The next portion of the analysis 
addresses the upwash spectrum and develops simplifications 
that lead to a tractable, computer-implementable form of 
Equations 22 and 23. 

The expression (Equation 23) just derived for the gener- 45 
alized upwash spectrum S ww (K,K',oo) is written in terms of 
vector wavenumber arguments. To proceed with its analysis, 
a coordinate system is chosen to represent the wavenumber 
vectors, and then the frozen turbulence assumption (Taylor’s 
hypothesis) is applied. Then, before re-inserting S ww , into 
Equation 22, the wavenumber coordinates are rotated into 50 
Glegg’s cascade system to facilitate manipulation of the 
noise equation. 

For reasons that will be apparent hereinafter, the preferred 
coordinate system for dealing with the turbulence inhomo- 
geneity is the x 0 , y 0 , z 0 system (FIG. 5), as illustrated in 55 
FIGS. 8(a) and 8(2?) . Turbulence is convected along the x 0 
axis. To ensure convergence of the Fourier integral, consider 
initially only a block of the flow, as indicated in FIG. 8(a). 
The streamwise and time extents of the flow source region 
are connected by the convection speed U 0 . Eventually, the 
extent of the source region is allowed to approach infinity, 60 
and the desired sound power is found by integrating the 
intensity over only the portion of the stator cascade corre- 
sponding to the stator vane of interest. In the z 0 direction, the 
source region is taken to be finite; the final noise formulas 
include an integral over this variable. 65 

To begin with, the transform of the velocity field given by 
Equation 15 is repeated below: 


Substitution back to the expression for S ww in Equation 23 
gives: 

2 JT , 1 (28) 

^ - k x U 0 )6(o> - *;t/o)yy 

ffffff myo ’ %) * ( *'’ y ' 0 ’ x ci ' M+Vo+^^^0, 

£ -t{ k x x ^yyo^z^dixdyodzodx'dyQdzo 


In Equation 28, the field points for the w’s are expressed 
independently. For extracting the covariance function below, 
the second field point is expressed based on a vector to the 
first plus a separation vector s via the following change of 
variables: 

X'=x+s x , y’ 0 =y 0 +s y , z' 0 =z 0 +s z (29) 


to get: 


= <5(oj - k x U 0 )S(oj - k' x U 0 ) x 

ffffff- 


(30) 


(2 tt) 5 T 

w(x, y 0 , zo)w{x + s x , yo + s y , zo + s z ) X 

* [( k x~ k x )x+{ky-k’ y ) y 0 +(k z -k' z )z 0 ) ] 

e~^ k x Sx+k y s y +k z s z)dxdy 0 dzods x ds y ds z 


Given the arguments of the two delta functions, then set 
k ;c =k , Ar , rearrange, and take the expected value (stochastic 
average) of both sides to get: 
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(S ww ) = S(a> -k x U 0 )S(co-k' x U 0 )x 


(31) 




_Uo_ 

(2tt) 5 

w(x + s x , >1 0 + Sy, Zo + s z ))dx x 

e -i(k^s x +k' y s y +k' z s z ) e i(ky -k' y ) y 0 £ i{k z -k' z ]ZQ 

d yodzods x dSy ds z 


10 


R(zo 


’ s)= iff { 


{w(x, ;y 0 , Zo)w(x + s x , y 0 + Sy, Zo + s z ))dx 


U 0 

(S ww ) = ~-6(o) - k x U 0 )6(oj - k’ x U 0 )6(k y - k' y ) x 


15 


There now follows a critical argument regarding the quantity 
in square brackets [0] in Equation 31. Note that <w(x,y 0 ,z 0 ) 
w(x+s JC ,y 0 +s 3 ; ,z 0 +sJ> is the covariance of the upwash veloc- 
ity component. At this point, it is totally general; that is, a 
function of the position vector (x,y 0 ,z 0 ) and the separation 
vector (s^s ,s*). Because U 0 T is the extent of the block of 
flow in the x direction (FIG. 8 (a)), the quantity in square 20 
brackets [0] in Equation 31 is the streamwise average of this 
covariance. If T is chosen to give an integral number of wake 
passages, and if block of flow is extended in the y 0 direction, 
then this average over x will be independent of y 0 . With the 
average covariance now independent of both x and y 0 , it can 
be considered to represent the average at the stator face (still 
a function of radius z 0 ). This function is named R according 
to: 


(32) 


(33) 40 


14 


so that: 


(Sw) = UqS((x) - k x U 0 )8{co - k' x U 0 )8(ky - k' y ) (35) 

2J K'y^o^dzo 


Before substituting this back to the acoustic equation, the k x , 
k y , k z wavenumbers must be expressed in the cascade 
system. This is accomplished by writing K-x in the cascade 
system and also in the x 0 , y 0 , z 0 system and equating them 
(because K-x is a scalar and therefore invariant under 
coordinate transformation): 


k^o+k^o+k^ff^oX+ay+vz 


(36) 


25 


We can write x, y, z in terms of x 0 , y 0 , z 0 , substitute on the 
right of Equation 36, and collect terms as follows: 

k x x 0 + k y y 0 + k z zo = (yocosi// s + vsim// s )x 0 + (37) 

+ arcos^/ + vsin^/cos^Jyo + 
(-yocosif/isiiuf/ s + arsing + vcos^/cos^ s )£o 


30 


However, despite the integration variable, the average could 35 
be taken across the stator face at constant z 0 . With the 
covariance function independent of y 0 , the y 0 integral yields 
a delta function and Equation 31 becomes: 


Equating coefficients of the space variables leads to the 
desired wavenumber transformation: 

^Jt = Yo C0S ‘ l l , s +V s i n ^5 

A^-Yq sin 'll?/ sin %+a cos i|v+v sin cos 

k=- Yq cos sin %-a sin cos cos (38) 


With these transformations, the desired goal of changing the 
represention of a band of turbulence oriented at constant 
radius in the duct is achieved. With the wavenumbers of 
Equation 38 substituted into Equation 35, S ww is now 
represented in the cascade coordinates where it runs diago- 
nally across the vane leading edges. This is the form required 
for manipulation as follows. 

Equation 22 can be written out with the wavenumber 
differentials shown explicitly as follows: 


pss e 




(39) 


55 


In the absence of the z 0 dependence, the quantity in square 
brackets would be the three-dimensional turbulence spec- 
trum for homogeneous turbulence (See Landahl, M. T. et al., 
“Turbulence and Random Processes in Fluid Mechanics”, 
Cambridge University Press, 1992). Hence, an analogous z 0 
dependent spectrum is defined according to: 


<p ww (Zo, r ) = ^3 fff R (Z0’ 


60 


(34) 


65 


This is the intensity spectrum as a function of the x, y, z 
space variables. Acoustic power is its integral over y and z 
as follows. For z, integrate over all z, recalling that the 
source is finite in that direction (because, temporarily, it is 
finite in the y direction). For the y direction, integrate over 
a distance corresponding to the number of vanes V in the 
cascade. The derivation is simpler if the z integration is 
carried out first, then inserting the expression for S H VV , and 
then performing the y integral. The z integral simply yields 
a delta function so that: 



15 
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p o n 2 c 4 0 J 

pss e 


ffffff—ZZ 


4^1 


(40) 


v') x 


10 


Then, when Equation 35 is inserted for with the 
wavenumbers transformed via Equation 38 and the v 1 inte- 
gral is performed, the result is: 


/■ 


I± dz =*!p^.x 

pSS e 


(41) 


ffffff- 


U 0 (y 0 cosijf s + vsim{/ s )d[a) - U 0 (y' 0 cosif/ S + vsim/r s )] X 


S[-(7o ~ 7o)si#,sin^ s + (a- T)cos^] x fi^Zo, K') 




k k > yiz-7k 

e^ x t- l '^-^ h ^- M) ^- M ' ) ' ylh Uy 0 dadvdy'Ja' 


30 


Note in the argument of the first delta function that U 0 cos 
i|^=U and U 0 sin %=W so that it becomes 6(co-y 0 U-vW). 
This is equivalent to 6 [y 0 -(co-vW)/U]/U enabling the y 0 and 
Yo integrals. As a result, in the following it is implicit that 
Y 0 =Y , o =c °g/U, where cog is Glegg’s shifted frequency co-vW. 
Equation 41 reduces to: 


/■ 


Itdz = ! 


(42) 


fisS e U 2 


Iff* - a')cosi//i] ^ (, 


<Pww(Zo, 


dzo X 


ZZ® 


£#ZT(A£)D(A£) 


k k’ > J K e~fk 

e m^>-^^^^m dadvda ^ 


40 


45 


50 

The remaining delta function can be re-written as 6(a-a')/ 
cos enabling the a’ integration: 



fisSgUcOSlf/iCOSlj/t 


Iff ' 


<pww(Zo, K')dZo x 


ZZ 


ffffiZT(A?)D(Aff) 

V *e - fk 


(43) 

55 


P>-y*i^w-*'w\ dadv 60 


where U=U 0 cos ty s . To complete the area integration, a 
coordinate r\ parallel to the leading edges in the z=0 plane 65 
is defined. Note that x-yd/h does not vary in this direction. 
Then, integrate over a distance equal to V vane gaps: 


f 


J2n(k-k' )- 


s drj = Vs6 kk 


(44) 


The k' summation can be completed via the Kroneker delta. 
With the primed wavenumbers and k' equal to their 
unprimed counterparts, the primed and unprimed Vs are 
equal and the exponential disappears completely. 

The sound power emanating from the V vanes is: 


n % = 


(45) 


Vp 0 7T 2 C 4 (j) 

fis e Ucosif/tcosif/ s 



E \£Dg t£ 


dadv 


This is the principal theoretical result of this analysis, and 
can be interpreted as follows, is the spectrum of sound 
power scattered upstream (+) and downstream (-) by the 
stator vane cascade as a result of interaction with turbulence 
described by the spectrum Thus, to compute sound 
power in any frequency band, integrate according to: 

n^ b = ir b nid(o (46) 

J(u a 


where the factor of 2 arises from the fact that the integration 
is only performed over positive frequencies. Equation 45 is 
the integral over the alpha a and v wavenumbers of the 
product of the turbulence spectrum (in square brackets) and 
the power response spectrum of the cascade. The turbulence 
wavenumber K has three components: y 0 , a, and v. The first 
of these, y 0 > is tied to the frequency by Yo=( w_v W)/U. The 
power response spectrum includes the factor in front of the 
integrals in Equation 45 and the k-summation. In perform- 
ing the wavenumber integral, the ranges of a and v are 
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limited to cut on waves; this amounts values of the wave- 
numbers leading to a positive argument of the square root in 
Equation 45. Explicit formulas for this range are given 
hereinafter. 

Equation 45 is not in the best form for computation 5 
because the response spectrum, which is time consuming to 
compute, is inside a double integral and a summation. This 
can be remedied by changing integration variables via the 
following procedure. First, move the summation outside the 
integrals in Equation 45 and make a slight change in notation 10 
by showing the components of the turbulence wavenumber 
explicitly as follows: 


18 

In Equation 50, the turbulence spectrum is expressed 
symbolically as dependent on the wavenumber components 
Yo, cx, v aligned with the vanes as in Glegg’s theory. Since 
lean and sweep are arbitrary, any turbulence spectrum used 
is likely to be based on other coordinates. For this 
application, LiepmamTs spectrum is used: 

_2^A 5 k\+k\ < 51 > 

” " 2 ([1 + A 2 (&i +&2 +^ 3>]) 3 


[£ = 


V p o 7t 2 C A OJ 
fis e UcOSl{/[COSlf/, 


[J~ <pww(Zo', Jo, CH, v)dZo\ 


V ' K e ~ fk 


dadv 


(47) 


There is an a integral for each k in the sum. The origin of 
each of these integrals can be shifted by noting in Glegg’s 
formulas that a only appears in the interblade phase angle 
a=y 0 d+ah. Furthermore, a appears only in the combination 
a-2 Jtk. Thus, if the integration variables change via: 

a-2 %klh=a' 

o-2 Jtk in the cascade response formulas changes to Yod 
+cx f h. With these changes, Equation 47 becomes: 


This applies to the upwash component which is associated to 
the “3” direction and whose mean square value is w 2 . A is 
the integral scale of the turbulence. To relate the k 1 , k 2 , k 3 
components to Yo? v > use the same technique applied 
earlier where Kx was equated in two different coordinate 
systems. One coordinate system is Glegg’s x, y, z system. 
For the other, the system is used in which one coordinate is 
aligned with the flow and another coordinate is aligned with 
a vane normal. This is the primed system illustrated in FIG. 
5. Thus: 


(48) 


25 


II' = - 


fiSeUcOSlf/tCOSlj/, 


a ' + Zxk/b’ v ) d Zo] 


V ^ - fo 


(49) 


da'dv 


35 

Since the cascade response term no longer depends on the k kpc'+k^'+^z'^ox+ay+vz (52) 

index, the expression can be rewritten in the final form: 


Yli = 


1 3s e Ucosij/icosij/ s 


r r ig0 z[fi 

J J V^ _ /o k 


<Pww(z 0 \ Jo, ® + 2 nk / h, v)dio \dadv 


(50) 


where the prime on the dummy variable a' is dropped. 45 
Programmed in software in this form, the inner loop has a 
summation on the turbulence spectrum that is more efficient 
than recomputing the cascade function per Equation 45. 
Equation 50 is the most general form of the noise theory for 
inhomogeneous turbulence impinging on a swept and/or 50 
leaned stator vane. It is provided as the basis for a noise 
prediction system for broadband rotor/stator interaction. 

Use of the general noise prediction theory is illustrated 
herein by a simple application of the Liepmann turbulence 
spectrum for homogeneous turbulence. The Liepmann tur- 
bulence spectrum is one of the classic spectrum formulas for 55 
homogeneous/isotropic turbulence. See Liepmann, H. W., 
“Extension of the Statistical Approach to Buffeting and Gust 
Response of Wings of Finite Span”, Journal of the Aero- 
nautical Sciences, March 1955. This application illustrates 
how to integrate over a turbulence field, and how to trans- 60 
form wavenumbers for turbulence represented in different 
coordinate systems. This application serves as an illustration 
of the use of Equation 50. To apply Equation 50, initially the 
form of the turbulence spectrum must be addressed and then 
a non-dimensionalization is chosen suitable for program- 
ming and for plotting of sound power spectra in decibels. 65 
Finally, specific guidance is provided for the integration and 
summation ranges that are currently implied to be infinite. 


Note that on the left side the “1” direction is associated with 
x', the flow direction, and the “3” direction is associated with 
y', which is the coordinate normal to the vane surface. By 
default, the “2” direction is aligned with z'. By substituting 
transformation equations relating the unprimed and primed 
coordinates and equating coefficients of x', y', and z', the 
desired connection is made as follows: 

£i=Yo cos t|j 5 +v sin 

k 2 =-y 0 sin 4> s +v cos % 

k 3 = a (53) 

However, recall that Yo was ti e d to frequency via Yo = ( C0 ” 
vW)/U. By substituting this into the first of Equations 53 and 
applying the relations U=U 0 cos % and W=U 0 sin ty s , k 2 
reduces simply to: 

&i=cd/[/q 


(54) 
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The same expression for y 0 leads to: 


*2 = - 


—tan ifj s + - 

U 0 cos i)/ s 


20 

The radial integration is performed over the turbulence (the 
inside integral) in a strip sense so that that integral can be 
(55) approximated by T> 33 AR. The result for one strip is: 


Also note for use in Equation 53 that a+2 Jtk/h is required, 
rather than just a. Hence: 

k 3 -> a+2 nk/h (56) 

Now, the turbulence spectrum is non-dimensionalized based 10 
on duct radius R, with the result: 


pwrP = 


2n 2 VM r (x) Aw w 2 
fis e cos& t cos 2 iJ/ s U 2 


lp 0 a 2 RARa\j' J* F ± {a, v)^ <b 33 dadv 


(63) 


(^ 33 — W R ^33 


where 



(57) 

The quantity in square brackets has the dimensions of 
15 power; the remainder of the expression is dimensionless. 

(58) Power level in decibels is given by 


20 


and K 1 =k 1 R, etc. The non-dimensional wavenumbers can be 
written in the required form for use with <D 33 in Equation 58: 


PWIj 1 



PWR t 
PWRref . 


(64) 


K[ = WjM r 

v 

K 2 = - (w / M r )tamf/ S + — 

CCOSl// s 



- + — 

C CSCOS^ 


(59) 


25 where the power reference is 10 12 watts or 


30 PWR^ f =0.73756x 10- 12 ft lb/sec 


(65) 


The overbars imply non-dimensionalization as follows 

co=coi?/a v=vc ct=ac c=c/R s=s/c (60) 

and M, =U 0 /a is the Mach number associated with the total The resu ll i s a working formula that can be coded in 
velocity vector. software for this particular application: 


PWR^ 

PWRref 


2tt 2 V M r l 2 (DA(D 
J3s e cosiJ/icos 2 if/ s 


\AP 0 RARa 

PWRref 




<f> 33 dadv 


(66) 


To obtain power in any band of frequency, integrate the 
spectral density over co according to Equation 46. For small 
enough bandwidth, however, simply multiply by 2xAoo. The 
factor of 2 is required to account for negative frequencies. 
This integrated quantity is referred to as PWR, which, with 
the above definitions applied to Equation 50, becomes: 


where p 0 a 2 =1.4P 0 is utilized, 1.4 is the ratio of specific heats 
for air, and P 0 is the ambient pressure. 

l 2 = ^IUl 

is the normalized mean square upwash velocity. 


45 
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PWRt = 


(61) 


27t 2 V p 0 a 2 M r Ra w 2 
/j^COS^/COS 2 ^ U 2 


ojAoj 


J F ± (a, v)£ [J" <J>33 dzo 


\dadv 


where the cascade power response function (in non- 
dimensional form and without the constants leading Equa- 
tion 50) is: 


The remaining items to be considered for programming in 
software for this particular application are the limits on the 
integrals and sum in Equation 66. These are treated as 
follows. The double integral in Equation 66 is over the a,v 
55 wavenumber plane. There is an ellipse in that plane such that 
wavenumbers inside it correspond to cut on waves and 
wavenumbers outside correspond to cut off waves. Cuton is 
governed by the square root in Equation 62. When the 
6Q argument of the square root 


_ I<£cD(A£)| 2 

F~(a, v) = 

yj c 2 (k 2 - /p 2 ) 


(62) iE=k e 2 -/* 2 (67) 

65 is positive, the waves are cut on. The cuton boundary is 
given by E=0, and k=0, in which case Equations 5 and 6 
yield: 
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, v 2 (cr + KMd) 2 (68) 

K ~w — =0 


This is the ellipse in the cr,v wavenumber plane shown in the 
sketch below. 
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This can be re-arranged and non-dimensionalized for pro- 
gramming as follows: 


ai = -(k / M) tanx - fP-K 2 - v 2 j co$x e 

a 2 = -(k I M)tanx + yj fl 2 * 2 - v 2 j cosx e 


v 


10 



Equation 68 can be solved for a as follows: 

cr = -KMd ± s e yj K 2 -v 2 j ft 2 

where, from Equations 3 and 5: 

w-vW 
K ~~^~ 

We will use the outside loop is used for v and integrated over 
the range of a corresponding to a 1 <a<a 2 where: 

cr l = -KMd -s e ^ k 2 -v 2 jfl 2 (71) 

cr 2 = -KMd + s e ^K 2 - v 2 / ft 2 ^ 

as shown in the sketch. The extremes of v occur when the 
square roots above are equal to 0 ; that is, for: 

45 

o)/a (72) 

Vl ’ 2 ~ Wfa + fi 


(69) 

30 

(70) 

35 


or in non-dimensional form for program coding in software 50 
as: 

(73) 

55 


where co=coc/a. Finally, since a=y 0 d+ah, the range of a 
corresponding to Equation 71 is given by: 

aci = -(■ y 0 +KM )^ - ^V/f 2 - v 2 //? 2 
or 2 = -{Jo+kM)^ + k 2 - v 2 / /? 2 


60 

(74) 

65 


-O) 

ft - M r smnf/ S 

O) 

ft + M r simf/ S 


where k = kc. 


In the software program code, the sum is initialized by the 
k=0 term. Then, a k loop runs from one to a large number, 
for example, thirty. The positive and negative terms are 
added in pairs until a convergence criterion is satisfied. 

In this exemplary application, Equations 58, 59, 62, 66, 73 
and 75 are coded in Fortran and executed by the computer 
or workstation 70. An exemplary input data set for this 
computer-implemented application of the noise prediction 
theory is illustrated in FIG. 9. The following inputs are 
selected by the user prior to computer program execution. In 
FIG. 9, line 1 is a comment line. In line 2, “G” is the 
gap/chord ratio, “theta” is stagger angle in degrees 
(measured from fan axis), “PsilD” and “PsisD” are the 
desired values for the lean and sweep angles, respectively, in 
degrees (in the “hat” convention of FIG. 6 ), and “Vanes” is 
the number of vanes in the stator cascade. In line 3, “Mr” is 
the stator inflow Mach number, “L” is turbulence length 
scale normalized by duct radius, and “I” is the rms upwash 
velocity normalized by the mean inflow velocity. Since the 
program code uses a homogeneous turbulence model in this 
exemplary application, “L” and “I” must represent values 
over a band of turbulence where the geometry, mean flow, 
and pitch averaged turbulence can be taken as uniform. 

In line 4, “Po” and “Ao” are the ambient pressure and 
speed of sound, “R” is the radius to the center of the band 
of turbulence, and “DelR” is its width. The units of “Po”, 
“R”, “DelR”, and “Ao” are foot pounds/sec. In line 5, 
constant bandwidth spectra or V 3 OB spectra can be com- 
puted. In line 6 , see the notes in FIG. 8 for specifying 
spectrum frequencies and bandwidths. Finally, in line 7, 
“DeltaAlf” and “DeltaNu” are step sizes for the a and v 
integrations. The values used (0.80) have been tested (by 
using smaller step sizes) and found adequate for this case. 
Other conditions may require smaller steps. 

In this exemplary embodiment, the computer program 
executes Equations 58, 59, 62, 66 , 73 and 75, together with 
the input file of FIG. 9, to compute the output sound power 
spectra (noise) emanating from the fan section using the 
user-input desired values for vane lean and/or sweep. Output 
sound power spectra are displayed on a visual display 
screen, and may also be written to a file. For example, a plot 
of sound power spectra may be created by reading the file 
into a Microsoft Excel spreadsheet and using its plotting the 
“frozen gust” assumption, or Taylor’s hypothesis. Turbu- 
lence at the stator face can be inhomogeneous and aniso- 
tropic. The source property required for a noise calculation 
is the radial integral of the pitch- wise average of the turbu- 
lence spectrum (or turbulence covariance function). The 
analysis contained herein gives the formulas for these opera- 
tions on the turbulence spectrum or covariance function. 

A special application of the theory has been developed by 
use of Liepmann’s turbulence spectrum (which is normally 
used for homogeneous/isotropic turbulence). This applica- 
tion has been described and illustrated herein with explana- 
tions of the required wavenumber transformations and 
details needed for software program coding. Equations are 
shown in the non-dimensional form used by the software 
program code. 

Although the present invention has been shown and 
described herein with respect to detailed embodiments 
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thereof, it will be understood by those skilled in the art that 
various other changes in the form and detail thereof may be 
made without departing from the broadest scope of the 
claimed invention in light of the teachings herein. 

Having thus described the invention, what is claimed is: 

1. A method in a computer system for determining an 
amount of noise at an inlet to a stator vane cascade in a gas 
turbine engine, the cascade including a plurality of stator 
vanes, the cascade being disposed downstream of a rotor 
blade array in the gas turbine engine, the method comprising 
the steps of: 

inputting a value indicative of an angular physical posi- 
tioning of at least one vane in the cascade; 

inputting a value indicative of an amount of turbulence 
emanating from the rotor blade array; and 

determining an amount of noise at the inlet to the stator 
vane cascade as a function of the value indicative of the 
angular physical positioning of the at least one vane in 
the cascade and as a function of the value for the 
amount of turbulence emanating from the rotor blade 
array. 

2. The method of claim 1, wherein the value indicative of 
the angular physical positioning of at least one vane in the 
cascade is indicative of a sweep angle of the at least one vane 
in the cascade. 

3. The method of claim 1, wherein the value indicative of 
the angular physical positioning of at least one vane in the 
cascade is indicative of a lean angle of the at least one vane 
in the cascade. 

4. The method of claim 1, wherein the value indicative of 
the angular physical positioning of at least one vane in the 
cascade is indicative of both a lean angle and a sweep angle 
of the at least one vane in the cascade. 

5. A computer system for determining an amount of noise 
at an inlet to a stator vane cascade in a gas turbine engine, 
the cascade including a plurality of stator vanes, the cascade 
being disposed downstream of a rotor blade array in the gas 
turbine engine, the computer system comprising: 

means for inputting a value indicative of an angular 
physical positioning of at least one vane in the cascade; 

means for inputting a value indicative of an amount of 
turbulence emanating from the rotor blade array; and 

means for determining an amount of noise at the inlet to 
the stator vane cascade as a function of the value 
indicative of the angular physical positioning of the at 
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least one vane in the cascade as a function of the value 
for the amount of turbulence emanating from the rotor 
blade array. 

6. The computer system of claim 5, wherein the value 

5 indicative of the angular physical positioning of at least one 

vane in the cascade is indicative of a sweep angle of the at 
least one vane in the cascade. 

7. The computer system of claim 5, wherein the value 
indicative of the angular physical positioning of at least one 
vane in the cascade is indicative of a lean angle of the at least 
one vane in the cascade. 

8. The computer system of claim 5, wherein the value 
indicative of the angular physical positioning of at least one 
vane in the cascade is indicative of both a lean angle and a 

15 sweep angle of the at least one vane in the cascade. 

9. A computer-readable medium containing instructions 
for determining an amount of noise at an inlet to a stator 
vane cascade in a gas turbine engine, the cascade including 
a plurality of stator vanes, the cascade being disposed 

20 downstream of a rotor blade array in the gas turbine engine, 
by: 

inputting a value indicative of an angular physical posi- 
tioning of at least one vane in the cascade; 

inputting a value indicative of an amount of turbulence 

25 emanating from the rotor blade array; and 

determining an amount of noise at the inlet to the stator 
vane cascade as a function of the value of the angular 
physical positioning of the at least one vane in the 

cascade and as a function of the value for the amount 
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of turbulence emanating from the rotor blade array. 

10. The computer-readable medium of claim 9, wherein 
the value indicative of the angular physical positioning of at 
least one vane in the cascade is indicative of a sweep angle 
of the at least one vane in the cascade. 

35 11. The computer-readable medium of claim 9, wherein 

the value indicative of the angular physical positioning of at 
least one vane in the cascade is indicative of a lean angle of 
the at least one vane in the cascade. 

12. The computer-readable medium of claim 9, wherein 

40 the value indicative of the angular physical positioning of at 
least one vane in the cascade is indicative of both a lean 
angle and a sweep angle of the at least one vane in the 
cascade. 





